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Abstract. We investigate the allowed range of reheating temperature values in light of the 
Planck 2015 results and the recent joint analysis of Cosmic Microwave Background (CMB) 
data from the BICEP2/Keck Array and Planck experiments, using monomial and binomial 
inflationary potentials. While the well studied cjP' inflationary potential is no longer favored 
by current CMB data, as well as with p > 2, a potential and canonical reheating 
{wre = 0) provide a good fit to the CMB measurements. In this last case, we find that the 
Planck 2015 68% confidence limit upper bound on the spectral index, Ug, implies an upper 
bound on the reheating temperature of Tre ^ 6 x 10^*^ GeV, and excludes instantaneous 
reheating. The low reheating temperatures allowed by this model open the possiblity that 
dark matter could be produced during the reheating period instead of when the Universe is 
radiation dominated, which could lead to very different predictions for the relic density and 
momentum distribution of WIMPs, sterile neutrinos, and axions. We also study binomial 
inflationary potentials and show the effects of a small departure from a potential. We 
find that as a subdominant cfP' term in the potential increases, first instantaneous reheating 
becomes allowed, and then the lowest possible reheating temperature of T^e = 4MeV is 
excluded by the Planck 2015 68% confidence limit. 
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1 Introduction 

Dark matter candidates, such as Weakly Interacting Massive Particles (WIMPs), sterile 
neutrinos, and axions are produced before Big Bang Nucleosynthesis (BBN). In order to 
preserve the thermal history of the Universe, including BBN and all subsequent events, the 
lower limit imposed by BBN on the highest temperature of the radiation dominated epoch 
(in which BBN happens), i.e. the reheating temperature, is 4MeV [1] (see also [2]). Notice 
that, for example, the standard freeze-out temperature of thermal WIMPs is Tfo — mj20, 
thus WIMPs with mass m > 80MeV are produced at temperatures above the BBN limit. 

We do not know the history of the Universe at temperatures above the BBN lower 
limit. Thus, in order to compute the relic abundance of dark matter candidates, assumptions 
must be made about the history of the Universe prior to BBN. The standard assumption is 
that the Universe was radiation dominated up to very large temperatures. However, there 
are other non-standard pre-BBN cosmological models, such as low reheating temperature 
models, that make different assumptions about the history of the Universe. These models 
make predictions of the relic density (see e.g. Ref. [3]) and momentum distribution (see e.g. 
Ref. [4]) of dark matter candidates that can be very different from the predictions made 
under the assumption that the dark matter is produced when the Universe is dominated by 
radiation. 

Recent experiments studying the Cosmic Microwave Background (CMB), such as the 
Background Imaging of Cosmic Extragalactic Polarization 2 (BICEP2) and Planck satellite 
experiments, are searching for evidence of primordial gravitational waves generated during 
inflation. While the BICEP2 experiment claimed the potential detection of a high level 
of primordial tensor modes in the polarization pattern of the CMB [5], the expected level 
of tensor modes shifted downward as the contribution due to foreground dust was better 
constrained. Indeed, a recent joint analysis of BICEP2/Keck Array [6] and Planck data has 
put an upper limit on the tensor-to-scalar ratio of amplitudes of perturbations produced 
during inflation at the scale (so called pivot scale) k = 0.05 Mpc“^, [7] 

r < 0.12 (1.1) 

at the 95% confidence level. This limit can be combined with the already measured amplitude 
of the scalar perturbations to find a bound on the energy scale of the inflationary potential 
V of < 1.2 X 10^6 GeV. 
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Inflation is a period of rapid expansion of the Universe before it became dominated 
by radiation. The usual inflationary models consist of a scalar field slowly rolling down an 
almost flat potential (see e.g. [8, 9]). Reheating happens when the slow-roll period ends and 
the energy of the inflaton held is converted into radiation (see [10] for a review). If reheating 
is assumed to be instantaneous, then the reheating temperature, T^e, is given by 

^ ^ '^9.{Tre)Tt (1.2) 

Assuming that the degrees of freedom come from the particles in the Standard Model (g* ~ 

107 for T > 175 GeV), Eq. 1.2 gives T^e < 7 x 10^^ GeV. If reheating is instantaneous, the 
Universe would be dominated by radiation up to this very large temperature. 

However, the physics of reheating is highly uncertain, and the reheating period may 
be extended. During this extended period, the energy density of the Universe decreases due 
to the expansion of the Universe, and is finally converted to a radiation bath at a smaller 
temperature than that given by Eq. 1.2 (see e.g. [11-15]). This opens up the possibility that 
the reheating temperature may not be much larger than the limit of 4 MeV imposed by BBN. 

Here we explore the possible range of reheating temperature values implied by some 
simple inflationary and reheating models, in light of the Planck 2015 results and the recent 
BICEP2/Keck Array and Planck joint analysis. In particular, we want to test whether T^e 
can approach the BBN limit in any of these models. We first consider the archetypal models 
of inflation, large field inflation models, which are characterized by a monomial potential. 

The simplest model, V = [16], has received ample attention recently because it fitted 

the previous data well [14, 17]. Thus, we consider similar monomial (V oc (j/') inflationary 
potentials and binomial {V oc + b(p'^) inflationary potentials, where one of the two terms 
is dominant. We consider inflationary models with potentials of Xcj)^ (motivated by 

axion monodromy models [18]), and (as in chaotic inflation [9]), and for the first time, 
as a modification of the monomial potential, cf)^ + for general p and q, but b(j)^~^ <C 1. 

The reheating period can be described with an effective equation of state for a cosmic 
fluid Wre = P/p (see e.g. Ref. [11]), where P is the pressure of the fluid, p is the energy 
density, and Wre is the equation of state parameter. The conservation of energy equation, 
d(/?a^) = —Pda^, where a is the scale factor, then implies that the energy density scales as 
p oc Qr-3(i+J«re)^ ^ Unlverse dominated by the oscillations of a scalar field around the 

minimum of a monomial potential, (jP^ during reheating (which can be different than the 
potential during inflation), Wre = (p ~ 2)/(p + 2) [19]. In the canonical model of reheating 
[20-22] , the inflaton field oscillates about the minimum of a quadratic potential, and decays 
to relativistic particles. In this case, Wre = 0. As in Refs. [11-14], we use a range of Wre 
values to parametrize the physics of the reheating period. Inflation ends when Wre — ~l/3, 
so during reheating, the equation of state parameter must be greater than Wre ^ —1/3. In 
addition, for scalar fields p > |P|, thus Wre < 1- In this paper, we consider models with 
Wre = 0, Wre = 1/3 (corresponding to U oc during reheating), and an intermediary value, 

Wre = 1/6, as an approximate sampling of the reasonable range of Wre- For example. Ref. [23] 
find values of Wre between 0 and 0.25, including preheating. We also consider two extreme 
exotic models with Wre = ~l/3 (as at the end of inflation) and Wre = 2/3 (corresponding to 
V oc 

The characteristics of inflation are already constrained by measurements of CMB anisotropies 
and large scale structures originating from quantum fluctuations of the inflaton and gravi¬ 
tational fields during inflation(see e.g. [24]). In addition to the BIGEP2/Keck Array and 
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Planck joint analysis bound on r, given in Eq. 1.1, we use the Planck 2015 measurement [25] 
of the spectral index 

n* = 0.9655 ± 0.0062, (1.3) 

and the primordial scalar amplitude 

ln(10^°^s) = (3.089 ± 0.036), (1.4) 


to constrain the shape of the potential. The measurements of and Ag are given here with 
68% confidence limits. During the slow roll phase, n^, Ag, and r can be related to the height 
of the inflationary potential, V, and its first and second derivatives (V^ and V^^) by 


- 1 ~ 


Ml 

P2 




(1.5) 


and 


~ 127r2M|’V| ’ 

v2 

(1.7) 


where Mp = 2.4 x 
parameters ey and 


10^^ GeV is the reduced Planck mass. 
W, given by 



During the slow roll phase, the 


( 1 . 8 ) 


and 


T]v = Mp 


V 


(1.9) 


are much smaller than 1. 

We will compute Ug and r in terms of the reheating temperature for a given inflationary 
potential and reheating model (parametrized by Wre)- In order to do so, we start by relating 
Tre with Nk and Pend, where Nk is the number of e-folds of expansion of the scale factor of 
the Universe between the moment when the pivot scale, k, exits the horizon and the end 
of inflation, and pend is the energy density at the end of inflation. This relation depends 
on the reheating model (namely Wre)- Assuming that the energy density does not change 
appreciably during inflation, pend is given by the product rAg. We then compute Nk in 
terms of Ug and r for each particular inflationary potential. Our work extends the analyses 
carried out in Refs. [11-14, 17, 26] by including, for the first time, an analytical treatment 
of binomial inflationary potentials, as well as monomial potentials, by taking a larger range 
of values of Wre, and by including the BICEP2/Keck Array and Planck joint analysis bound 
on r as a constraint. Our approach is complementary to that of Ref. [27], which studies the 
relationship between the inflaton decay rate and Tre- 

We find particularly interesting that, while the well studied cfP' inflationary potential re¬ 
produces the value of Ug favored by Planck, the upper limit on r found by the BICEP2/Keck 
Array and Planck joint analysis cannot be satisfied with this potential by any of the re¬ 
heating models we consider, in agreement with Fig. 21 of Ref. [25]. The quadratic potential, 
therefore, is no longer favored; a potential, however, does provide a good fit to the CMB 
measurements. For this potential, the Planck 2015 68% confidence limit on Ug and the upper 
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bound on r imply an upper bound T^e ^ 6 x 10^*^ GeV for the reheating temperature, assum¬ 
ing canonical reheating. Only exotic reheating models with Wre < 0, such as Wre = ~l/3, 
allow for higher reheating temperatures, with a inflationary potential. If the Universe 
has a low reheating temperature, as allowed by this model, the possibility is opened that 
dark matter could be produced during the reheating period instead of when the Universe is 
radiation dominated, which could lead to very different predictions for the relic density and 
even momentum distribution of WIMPs, sterile neutrinos, and axions from the predictions 
assuming that dark matter is produced when the Universe is radiation dominated. 


2 Relationship between Nk and Tre 


In this section, we rederive a relation between the number of e-folds, Nf^ = ln(aend/®fc)) 
where Uend and are the scale factors at the end of inflation and when the pivot scale exited 
the horizon, respectively, and the reheating temperature, T^e in terms of known quantities k, 
r, As, Oeq, and peq (the last two are the scale factor and energy density at matter-radiation 
equality, when the energy densities of matter and radiation are equal). We do this by tracing 
the dilution of the energy density of the Universe during the reheating phase. 

During reheating, the energy density, p, scales as p oc a-K^+'^re) ^ where Wre is the 
equation of state parameter. This means that 

Pre \ ^end J 


where Pend^ Pre and a^nd: o,re are the energy density and scale factor at the end of inflation and 
at the end of reheating, respectively. After reheating comes an epoch of radiation domination 
(with p oc a~^), which lasts until the Universe becomes matter dominated at matter-radiation 
equality. Then 

^=f^\ . ( 2 . 2 ) 

Peq V Ore / 

Combining Eqs. 2.1 and 2.2, we find 


1 — 1 

aend _ ( Peq\ / Peq^. I2(“re+l) 

\ Pend J \ Pre J 


(2.3) 


Since aend = afcexp(Aifc), we can replace Uend in Eq. 2.3 and solve for Nk, 


Nk=ln(^] + —^ - In ( \ In (2.4) 

V Ofc y 3(1 -h Wre) \PendJ l2{Wre + 1) \PreJ 

This expression is valid for any inflationary potential. This potential is needed to relate pend 
to the energy density at the pivot scale pk- An equation similar to Eq. 2.4 is given originally 
in Eq. 15 of Ref. [12], and again (but solving for We = ^^{cLre/o-end) and applied only to the 
case of a monomial inflationary potential) in Eq. 11 of Ref. [14]. 

The energy density at matter-radiation equality is peg ^ ‘2'Prad = 2 ((7r^/30)p*^egE4) , 
where Teq = {ao/aeq)To and = 3.36 are the temperature and number of relativistic 
degrees of freedom at matter-radiation equality, respectively. The energy density at the end 
of the reheating epoch is related to the reheating temperature by pre = (7r^/30)p*^reEye' 
The number of relativistic degrees of freedom at the end of reheating, p*,re depends on 
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Figure 1. for pivot scale k = 0.05 Mpc“^ as a function of the reheating temperature T^e for 
Wre — 0 (black), Wre = 1/6 (red), Wre = 1/3 (blue), Wre = —1/3 (purple), and Wre = 2/3 (green), 
using Us = 0.9655, Ag = 2.2 x 10“®, and a^qja^ = 3360. To make this figure, we assumed that the 
inflationary potential is quadratic. 


the reheating temperature, and is approximately given by = 107, 90, and 11 for 
Tre > 175 GeV, 175 GeV > T^g > 200 MeV, and 200 MeV > T^g > IMeV, respectively, 
where T ~ 200 MeV roughly corresponds to the temperature of the QCD phase transition, 
and T < 175 GeV is when the temperature drops below the top quark mass [28]. 

There are still two unknown parameters, and Pend^ in Eq- 2.4. can be related 
to k, r, and As because at the moment the comoving scale k exits the horizon, the Hubble 
expansion rate, Hk cs ^[tt"^/ 2)MprAs, is related to k and ak by Hk = k/au- Although we 
do not know Pend-, we do know the energy density when the pivot scale exits the horizon, 
Pk = 3MpH‘^ ~ /2)MprAs. Given a particular inflationary potential, we can calculate 
the relationship between r and n^, which will allow us to find the values of Hk and pk by 
using the Planck measurement of Ug given in Eq. 1.3. If the energy density during inflation is 
approximately constant, then pg^^ ~ pk- As we will show later, this approximation introduces 
an error of no more than 6% to the value of Nk in the models with Wre = 0 that we consider 
in subsequent sections. 

Fig. 1 shows Nk as given in Eq. 2.4 as a function of the reheating temperature T^g for 
Wre = 0 (black), Wre = 1/6 (red), Wre = 1/3 (blue), Wre = —1/3 (purple), and Wre = 2/3 
(green), for k = 0.05Mpc”^, the central values of the measurements of n* and Ag, given 
in Eqs. 1.3 and 1.4, and aeqjao = 3360 [28]. To make this hgure, we assumed that the 
inflationary potential is quadratic. We note that taking values of Ug and Ag at the edge of 
the 68% confidence limit ranges given in Eqs. 1.3 and 1.4, respectively, will only change the 
calculated value of Nk by ANk ^ 0.1. The plot ends at T^g ~ 7 x 10^® GeV because this 
reheating temperature occurs when the inflaton field decays instantaneously into radiation, 
thus Pend does not redshift before being transformed into radiation. The case of Wre = 1/3 
gives a straight line because the expansion rate in the reheating epoch is the same as in 
a radiation dominated epoch. In this case, it does not matter when the Universe becomes 
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radiation dominated, because it expands in the same way, i.e. p oc a~^, both before and after 
the transition. 

Fig. 1 makes it clear that if the reheating temperature is known, the number of e- 
folds from when the pivot point exits the horizon to the end of inflation can be calculated, 
given an inflationary potential and a reheating model. As previously mentioned, the lowest 
possible reheating temperature is Tre — 4MeV. As can be seen in Fig. 1, this means that 
Nk ~ 42.5 is the lowest allowed number of e-folds between the time when the pivot scale of 
k = 0.05 Mpc“^ exits the horizon and the end of inflation for Wre = 0, assuming a quadratic 
inflationary potential. Lower bounds on Nk for other values of Wre can be found similarly. 
The upper bound for all reasonable reheating models with Wre < 1/3 is A'fc < 57. 

We made Fig. 1 under the assumption that the energy density is approximately constant 
during inflation. If the energy density does change during inflation, i.e. pk > Pend-, we must 
add the term 

1 ( 2 -^) 
3(1 + Wre) \Pend J 

to the right side of Eq. 2.4. In the models we consider in the subsequent sections, the ratio of 
the energy density when the pivot scale leaves the horizon to the energy density at the end 
of inflation varies from Pk/Pend — 15 to Pk!Pend — 3700, which leads to a shift of ^Nk — 0.9 
to AA^fc ~ 2.7 for Wre = 0. A similar estimation can be made for other values of Wre^ and 
AAfc is always small. 


3 Constraints on Tre from Us and r 


In this section we find expressions for the spectral index, n*, and the tensor-to-scalar ratio, 
r, for monomial {V oc (fT) and binomial (V = 4^ + bcj)^) potentials, as a function of Nk- We 
then use the ranges of r and n*, given in Eqs. I.l and 1.3 respectively, to derive constraints 
on Nk, and subsequently Tre- 


3.1 Monomial potentials 


The simplest models of the inflaton potential are monomial potentials. These potentials have 
the form 

(jtY 


V{(j)) = XMp 

The slow roll parameters ey and rjv ai'c given by 


V Mp J 


Mj, 


_P^ (Mp 

2 V 


and 




pv = Ml^ =p{p-i) 


(f) 


(3.1) 


(3.2) 


(3.3) 


At the end of inflation, i.e. when the slow roll approximations break down, the first slow roll 
parameter is generally taken to be ey = ^end — 1- Thus the value of the inflaton field at the 
end of inflation is 

(j^end P 


Mp c/2eend 


(3.4) 
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When the slow roll approximation holds, 


Nk ~ 


I [■4>k 

J<P<,ud 



(3.5) 


where Mp is the reduced Planck mass, V is the inflaton potential, is the derivative of the 
potential, and 4>k and (pend are the values of the inflaton field when the pivot scale exits the 
horizon and at the end of inflation, respectively. 

Using Eq. 3.5, we can then find an expression for in the case of monomial potentials 


Nk 


1 

( 4>ky 

p2 

2p 

\MpJ 

‘^^end 


(3.6) 


This relation can be inverted to solve for the value of the inflaton field when the pivot scale 
exits the horizon 

and finally we find expressions for ey and ijv as a function of Nk 

ev = 

and 

W = 


P^end 

^^k^end P 

(3.8) 

2{p l')6end 
^Nkeend+P' 

(3.9) 


The spectral index, Us, and the tensor-to-scalar ratio, r, are given in terms of the slow 


roll parameters. 


ns — 1 ~ 2riv — 6ey, 

(3.10) 

and 


r ~ 16e\/, 

(3.11) 

which for monomial potentials gives 


— 2(j? -|- 2')€end 
~ ^Nkeend+P ’ 

(3.12) 

and 

16p6end 

rp 

^Nk€end T P 

(3.13) 


We can see from Eqs. 3.12 and 3.13 that for monomial inflationary potentials, the tensor-to- 
scalar can be related to the spectral index by 


n., - 1 = 


-{p + 2)r 
8p 


(3.14) 


Fig. 2 shows the spectral index, Ug (left vertical axis) and the tensor-to-scalar ratio, r 
(right vertical axis), measured at pivot point k = 0.05 Mpc“^ as a function of Tre for p = 2, 
and Wre = 0 (black), Wre = 1/6 (red), Wre = 1/3 (blue), Wre = —1/3 (purple), and Wre = 2/3 
(green). To find Nk, we used the same values used in Fig. 1. The light red area marks the 
68% confidence limit range of the Planck 2015 measurement of Ug, while the light blue area 
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Figure 2. The spectral index, Ug (left vertical axis) and the tensor-to-scalar ratio, r (right vertical 
axis), measured at pivot point k = 0.05 Mpc“^ as a function of T^e for p = 2, and Wre = 0 (black), 
Wre = 1/6 (red), Wre = 1/3 (blue), Wre = —1/3 (purple), and Wre = 2/3 (green). To find Nk, we used 
the same values used in Fig. 1. The light red area marks the 68% confidence limit range of the Planck 
2015 measurement of Ug, while the light blue area marks the range of r allowed by the BICEP2/Keck 
Array and Planck joint analysis. The light purple area is where these two regions overlap. 


marks the range of r allowed by the BICEP2/Keck Array and Planck joint analysis. The 
light purple area is where these two regions overlap. 

We clarify that we used the central value of Ug = 0.9655 to calculate the value of the 
inflaton potential in order to find N^- is in turn used to calculate r and Ug. However, if 
Ug is allowed to vary between 0.87 and 0.97, as predicted by the models we consider in Fig. 2, 
the value of Nk is changed by no more than 0.5 (i.e. no more than 1.5%) for Wre = 0. 

We see that for the well studied (fP' inflationary potential, while the models reproduce 
the value of Ug measured by Planck, the upper limit on r found by the BICEP2/Keck Array 
and Plank joint analysis cannot be satisfied by any of these reheating models. Even for the 
exotic Wre = 2/3 model, the predicted value of r does not fall below r < 0.12 until the 
predicted reheating temperature drops below the BBN limit of 4 MeV. Thus we see that the 
(jp inflaton potential is no longer favored by current CMB data, in agreement with Fig. 21 of 
Ref. [25]. 

Figs. 3 and 4 are the same as Fig. 2, but for p = 1 and p = 4, respectively. We find 
that a (jP potential fits the CMB measurements better than the (jp model. We see that 
the 68% confidence limit on Ug implies an upper bound of T^e ^ 6 x 10^° GeV, assuming 
canonical reheating. Only exotic reheating models with Wre < 0, such as Wre = ~l/3, allow 
for higher reheating temperatures. In fact, with Wre = ~l/3, the 68% confidence limit range 
in Ug translates to the range 10^ GeV < Tre ^ 2 x 10^^ GeV for the reheating temperature. 
Instantaneous reheating is not allowed by the CMB measurements. Models with Wre = 1/6 
restrict the reheating temperature to Tre ^ 7x 10^ GeV, while models with Wre >1/3 predict 
values of ng that are outside of the Planck 2015 68% confidence limit range for all possible 
reheating temperatures. We find that for a potential, the predicted values of Ug and r are 
outside of the GMB limits for all possible reheating temperatures. 

In this section, we have presented bounds on Tre for models with power law index p = 1, 
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Figure 3. Same as Fig. 2, but for p = 1. 
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Figure 4. Same as Fig. 2, but for p = A. 


2, and 4, and equation of state parameter Wre = 0, 1/6, 1/3, —1/3, and 2/3. However, a 
similar analysis, making use of Eqs. 3.12 and 3.13, can be performed for models with other 
power law indexes and equation of state parameters. 

3.2 Binomial potentials 

We have seen in the previous section that instantaneous reheating is not allowed in the (j)^ 
model. We now examine analytically, for the first time, the impact of a small departure 
from this potential by studying binomial inflationary potentials. By adding a second term to 
the potential, we wish to find how much of a departure is necessary to allow instantaneous 
reheating, and then how much of a departure it takes before the lowest possible reheating 
temperature of 4 MeV is excluded. We thus consider models with inflationary potenials with 
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two terms, taking the form 


y((^) = \Mj 




^ \Mp J \Mp 

With this potential, the slow roll parameters are 

p + bq 


(3.15) 


_ 1 fMp\ 

n J 


Mp I 




q-p 


(3.16) 


and 


, Q-p 


W = 


(3.17) 


!'(?'-1)+ '>«(«-1) (jfe) 

The number of e-folds from when the pivot scale exits the horizon until the end of inflation 


IS 


1 r<t>k 


l + b(^ 


Mp I 


q-p 






(3.18) 


If the first term of the potential dominates over the second term during inflation, i.e. 

q-p 


1 > b 


(-)' 

\MpJ 


and assuming q/p ~ 0(1), Eq. 3.18 simplifies to 


nr 1 1 ^ 

^P 4end P ^P 

The solution to this integral is 
Nk ~ — 

q-p+ 2 


1 + 




q-p 



\Mp) 


4^end 

JPp 


dcj). 


q-p+2 


(3.19) 


(3.20) 


(3.21) 


The first term on the right side of Eq. 3.21 is the same as in the case of a monomial potential, 
and the second term is the first order correction. 

As in the case of the monomial potential, we want to find expressions for the slow roll 
parameters as functions of Nj-. If the first term in the potential dominates over the second, 
Eq. 3.16 simplifies to 


ev — 


2 



\Mp) 


(3.22) 


Rearranging this formula, and again using the approximation that the second term of the 
potential is subdominant, we find that 



2 





2ev 



\ / p2 \ (g-p+2)/2 

) \2wJ 


(3.23) 
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Using this expression in Eq. 3.21 and keeping only the first order correction, we find 


1 


P 


Nk--{^ + 2b ^ 


2p 2ev 


/^xU-p+2)/2 


\ 2 ey J 


--1 1 - 


1 


q-p + 2 


(3.24) 


A similar process can be followed for the second slow roll parameter, pv 




1 j pip — 1) , 1)V'^ P+2)/2 


VV 


+ b 


Vv 


qjq - 1 ) -pjp - 1 ) 

pip - 1 ) 


+ 


Hill 

q-p+2 


which is valid for all p 7 ^ 1. In the case oi p = 1, we find instead 

,^,2/(3-q) A n\3/(3-g) 

y VV J 3 




b nqiq-l) y 

vv J 


(3.25) 


(3.26) 


By inverting Eqs. 3.24 and 3.25 (or 3.26 for j? = 1) we can then find an equation for ey 
and r]v in terms of the parameter b and N^. To do this, we must use particular values of p 
and q. Once we solve these equations, we use Eqs. 3.10 and 3.11 to find the predicted values 
of Ug and r. 


3.3 Study of U oc (ji/Mp + 6(</>/Mp)^ potentials 

We now focus on the interesting case of a binomial potential with p = 1 and q = 2. We 
particularly wish to find out how much of a departure from a simple monomial potential 
is necessary to allow instantaneous reheating, and how much of a departure it takes before 
the lowest possible reheating temperature of 4MeV is excluded by the Planck 2015 68 % 
confidence limits. 

Before plotting Ug and r as functions of h is useful to know which values of b are 
allowed by the approximation that the second term in the potential is subdominant to the 
first. Eor p = 1 and q = 2, we need to satisfy 1 6 ((/>/Mp)^ during inflation. Keeping 

the leading term in Eq. 3.22, (c/i/Mp)^ ~ p^/(2ey) and using Eq. 3.24, even with at its 
maximum value of 57 (assuming canonical reheating), we find that b < 10“^ is enough to 
ensure that the second term in the potential is subdominant to the first. 

Eig. 5 shows the spectral index, Ug (left) and the tensor-to-scalar ratio, r (right), mea¬ 
sured at pivot scale k = 0.05 Mpc“^ as a function of T^e for a binomial potential with p = 1, 
q = 2, b = 10 “^, and Wre = 0 (black dashed), Wre = 1/6 (red dashed), Wre = 1/3 (blue 
dashed), Wre = —1/3 (purple dashed), and Wre = 2jS (green dashed). The solid black, red, 
blue, purple, and green lines show Ug and r for the corresponding monomial potential with 
p = 1. To find Nk, we used the same values used in Eig. 1. The red and blue shaded ar¬ 
eas mark the 68 % confidence limit region allowed by the Planck 2015 measurement of the 
spectral index and the BICEP2/Keck Array and Planck joint analysis upper bound on the 
tensor-to-scalar ratio, respectively. We are no longer able to plot both Ug and r on the same 
plot, since the relationship between the two now depends on the parameter b. As expected, 
when the second term in the potential is subdominant to the first, the dashed lines closely 
follow the solid lines for a monomial potential with p = 1, but are slightly offset in the direc¬ 
tion of the lines for a monomial potential with p = 2. As b increases, the dashed lines move 
increasingly farther from the solid lines, and toward the lines for a monomial potential with 

p = 2. 
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Figure 5. The spectral index, rig (left) and the tensor-to-scalar ratio, r (right), measured at pivot 
scale k — 0.05Mpc“3 as a function of T^e for a binomial potential with p = 1, q = 2, b = 10~^, 
and Wre = 0 (black dashed), Wre = 1/6 (red dashed), Wre = 1/3 (blue dashed), Wre = —1/3 (purple 
dashed), and Wre = 2/3 (green dashed). The solid black, red, blue, purple, and green lines show Ug 
and r for the corresponding monomial potential with p = \. To find Nk, we used the same values used 
in Fig. 1. The red and blue shaded areas mark the 68% confidence limit region allowed by the Planck 
2015 measurement of the spectral index and the BICEP2/Keck Array and Planck joint analysis upper 
bound on the tensor-to-scalar ratio, respectively. We are no longer able to plot both rig and r on the 
same plot, since the relationship between the two now depends on the parameter h. 


In particular, we want to understand how the limits on Tre change as h increases. The 
range of values of T^e with predicted values of r and Ug that fall within the ranges given in 
Eqs. 1.1 and 1.3, respectively, can be read off of Fig. 5. In the case of canonical reheating 
{wre = 0 ), with p = 1 , q = 2 , and b = 10 “^, the constraints are relaxed, compared to those 
of a monomial potential with p = 1 , because the addition of the second term in the potential 
lowers the predicted value of n^. In this case, we see that instantaneous reheating is allowed 
(as compared with an upper limit of T^e < 6 x 10^33 GeV for a pure monomial potential with 
p = 1). The lowest reheating temperature of 4Mev is also allowed in this case. 

Finally, we can find how much deviation from a monomial potential with p = 1 is 
necessary to allow instantaneous reheating, and how much of a depature it takes before 
the lowest possible reheating temperature is no longer allowed. We find that in the case 
of binomial potentials with p = 1 and g = 2 , as the subdominant term increases (i.e. h 
gets larger), instantaneous reheating is allowed for the first time when b = 5.3 x 10“^, when 
the predicted value of falls within the 68 % confidence limit range of the Planck 2015 
measurement. Assuming canonical reheating, the lowest reheating temperature of 4MeV 
is also allowed by this value of b. As the subdominant term continues to increase, the 
predicted value of n* at T^e = 4MeV first falls out of the 68 % confidence limit range of the 
Planck measurement when b = 1.4 x 10“2. This means that when b takes values between 
5.3 X 10“^ < b < 1.4 X 10“^, all possible values of the reheating temperature are allowed by 
the Planck 2015 68 % confidence limit range on Ug and the BICFP2/Keck Array and Planck 
joint analysis upper bound on r, assuming canonical reheating. 

While we have considered the interesting case of a binomial potential with p = 1 and 
q = 2, a similar analysis for inflaton potentials with different values of p and q can be carried 
out using the procedure described in Section 3.2. Corresponding limits on the reheating 
temperature of the Universe can then be found for these models as well. 
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4 Conclusions 


In this paper we rederived the relationship between the number of e-folds from when the pivot 
scale exits the horizon until the end of inflation and the reheating temperature. We find that 
in the case of canonical reheating, and a quadratic inflationary potential, the lower bound 
on Nk, measured at pivot scale k = 0.05Mpc“^, is > 42.5. In addition, for reheating 
scenarios with Wre < 1/3, IVfc cannot exceed 57. If our assumption that the energy density 
of the Universe remains constant during inflation does not hold, then the upper limit will be 
raised. For example, if the energy density falls by a factor of Pk/Pend — 3700 between when 
the pivot scale left the horizon and the end of inflation, as in the most quickly falling potential 
we considered, then the upper limit will be raised by ~ 2.7 e-folds, for Wre = 0. 

We then found the dependence of the spectral index, n^, and the tensor-to-scalar ratio, 
r, on Nk (and thereby Tre) for the archetypal models of inflation, large field inflation models, 
which are characterized by a monomial potential. We also studied, as a modification of the 
monomial potential, binomial inflationary potentials in which one term dominates over the 
other. We used the measurement of n* made by Planck 2015 (given in Eq. 1.3) and the upper 
limit on r given by the joint analysis of the BICEP2/Keck Array and Planck experiments 
(given in Eq. 1.1) to constrain the possible reheating temperatures of these models. We see 
that for the well studied (fP' inflationary potential, while the models reproduce the value of 
Us favored by Planck, the upper limit on r found by the BICEP2/Keck Array and Plank 
joint analysis cannot be satisfied by any of the reheating models we consider, in agreement 
with Fig. 21 of Ref [25]. Even for the exotic Wre = 2/3 model, the predicted value of r does 
not fall below r < 0.12 until the predicted reheating temperature drops below the BBN limit 
of 4MeV. Thus we find that the (jP inflaton potential is no longer favored by current CMB 
data. 

We see that a pP potential provides a good fit to the CMB measurements. We find 
that the Planck 2015 68 % confidence limit range on Ug implies an upper bound of Tre ^ 
6 X 10^^ GeV, assuming canonical reheating. Only exotic reheating models with Wre < 0, 
such as Wre = —1/3, allow for higher reheating temperatures. In fact, with Wre = —1/3, the 
68 % confidence limit range in Ug translates to the range 10^ GeV < Tre ^ 2 x 10^^ GeV for the 
reheating temperature. Instantaneous reheating is not allowed by the GMB measurements. 
Models with Wre = 1/6 restrict the reheating temperature to Tre x 10^ GeV, while models 
with Wre >1/3 predict values of Ug that are outside of the 68 % confidence limit range of the 
Planck measurement for all possible reheating temperatures. We see that for a 4>‘^ potential, 
the predicted values of Ug and r are outside of the GMB limits for all possible reheating 
temperatures. 

We then explored, for the first time, inflationary models with binomial potentials {V oc 
(jp + bcp'^) for general p and q, and use our results to find how much deviation from a monomial 
potential with p = 1 is necessary to allow instantaneous reheating, and how much of a 
depature it takes before the lowest possible reheating temperature is no longer allowed. We 
find that in the case of binomial potentials with p = 1 and q = 2, as the subdominant 

term increases (i.e. b gets larger), instantaneous reheating is allowed for the first time 
when b = 5.3 x 10“^, when the predicted value of Ug falls within the 68% confidence limit 
range of the Planck 2015 measurement. Assuming canonical reheating, the lowest reheating 
temperature of 4 MeV is also allowed by this value of b. As the subdominant term continues 
to increase, the predicted value of Ug at Tre = 4 MeV first falls out of the 68% confidence 
limit range of the Planck measurement when b = 1.4 x 10“^. This means that when b takes 
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values between 5.3 x 10“^ < b < 1.4 x 10“^, all possible values of the reheating temperature 
are allowed by the Planck 2015 68% confidence limit range on and the BICEP2/Keck 
Array and Planck upper bound on r, assuming canonical reheating. 

We find particularly interesting that for a monomial potential with p = 1, the 68% 
confidence limits on and the upper bound on r imply an upper bound T^e ^ 6 x 10^*^ for 
the reheating temperature, assuming canonical reheating. If the Universe has a low reheating 
temperature, like those allowed in this model, the possibility is opened that dark matter could 
be produced during the reheating period instead of when the Universe is radiation dominated, 
which could lead to very different predictions for the relic density and even momentum 
distribution of WIMPs, sterile neutrinos, and axions from the predictions assuming that 
dark matter is produced when the Universe is radiation dominated. 

Acknowledgements 

T.R. and G.G. were supported in part by the Department of Energy under Award Number 
DE-SC0009937. This research was also supported in part by the National Science Eoundation 
under Grant No. PHYl 1-25915 (through the Kavli Institute for Theoretical Physics, KITP, 
at the University of Galifornia, Santa Barbara, where G.G. carried out part of the work). 

References 

[ 1 ] S. Hannestad, What is the lowest possible reheating temperature?, Phys. Rev. D 70 , 043506 
(2004) [astro-ph/0403291]. 

[2] F. De Bernardis, L. Pagano and A. Melchiorri, New constraints on the reheating temperature of 
the universe after WMAP-5, Astropart. Phys. 30 , 192 (2008). M. Kawasaki, K. Kohri and 

N. Sugiyama, MeV scale reheating temperature and thermalization of neutrino background, 

Phys. Rev. D 62 , 023506 (2000) [astro-ph/0002127]. M. Kawasaki, K. Kohri and N. Sugiyama, 
Cosmological constraints on late time entropy production, Phys. Rev. Lett. 82 , 4168 (1999) 
[astro-ph/9811437]. 

[3] G.Giudice, E.Kolb, and A.Riotto Largest temperature of the radiation era and its cosmological 
implications, Phys. Rev. D 64 (2001) 023508 [arXiv:hep-ph/0005123]; N.Fornengo, A.Riotto, 
and S.Scopel Supersymmetric Dark Matter and the Reheating Temperature of the Universe, 
Phys. Rev. D 67 (2003) 023514 [arXiv:hep-ph/0208072]; C.Pallis Massive Particle Dacay and 
Cold Dark Matter Abundance, Astorpart. Phys. 21 (2004) 689 [arXiv:hep-ph/0402033]; 

G. Gelmini and P.Gondolo, Neutralino with the right cold dark matter abundance in (almost) 
any supersymmetric model, Phys. Rev. D 74 (2006) 023510 [arXiv:hep-ph/0602230]; M. Drees, 

H. Iminniyaz and M. Kakizaki, Abundance of cosmological relics in low-temperature scenarios, 
Phys. Rev. D 73 (2006) 123502 [arXiv:hep-ph/0603165]; G. Gelmini et al The effect of a late 
decaying scalar on the neutralino relic density, Phys. Rev. D 74 (2006)083514 
[hep-ph/0605016]. G. B. Gelmini, J. H. Huh and T. Rehagen, “Asymmetric dark matter 
annihilation as a test of non-standard cosmologies,” JCAP 1308 , 003 (2013) [arXiv: 1304.3679 
[hep-ph]]. L. Visinelli and P. Gondolo, “Axion cold dark matter in non-standard cosmologies,” 
Phys. Rev. D 81 , 063508 (2010) [arXiv:0912.0015 [astro-ph.GO]]. G. Gelmini, 

S. Palomares-Ruiz and S. Pascoli, “Low reheating temperature and the visible sterile neutrino,” 
Phys. Rev. Lett. 93 , 081302 (2004) [astro-ph/0403323]. G. Gelmini, E. Osoba, 

S. Palomares-Ruiz and S. Pascoli, “MeV sterile neutrinos in low reheating temperature 
cosmological scenarios,” JCAP 0810 , 029 (2008) [arXiv:0803.2735 [astro-ph]]. G. E. Yaguna, 
“Sterile neutrino production in models with low reheating temperatures,” JHEP 0706 , 002 
(2007) [arXiv:0706.0178 [hep-ph]]. T. Rehagen and G. B. Gelmini, “Effects of kination and 
scalar-tensor cosmologies on sterile neutrinos,” JCAP 1406 , 044 (2014) [arXiv:1402.0607 


- 14 - 



[hep-ph]]. L. Roszkowski, S. Trojanowski and K. Turzyski, “Neutralino and gravitino dark 
matter with low reheating temperature,” JHEP 1411, 146 (2014) [arXiv: 1406.0012 [hep-ph]]. 

[4] G. Gelmini and C. E. Yaguna Constraints on Minimal SUSY models with warm dark matter 
neutralinos, Phys. Lett. B 643 (2006) 241 [hep-ph/0607012]; G. B. Gelmini and P. Gondolo 
Ultra-cold WIMPs: relics of non-standard pre-BBN cosmologies, JCAP 0810 (2008) 002 
[astro-ph/0803.2349]. 

[5] P. A. R. Ade et al. [BIGEP2 Gollaboration] , “Detection of B-Mode Polarization at Degree 
Angular Scales by BICEP2,” Phys. Rev. Lett. 112, 241101 (2014) [arXiv: 1403.3985 
[astro-ph.GO]]. 

[6] P. A. R. Ade et al. [BIGEP2 and Keck Array Gollaborations] , “BICEP2 / Keck Array V: 
Measurements of B-mode Polarization at Degree Angular Scales and 150 GHz by the Keck 
Array,” arXiv:1502.00643 [astro-ph.GO]. 

[7] P. A. R. Ade et al. [BIGEP2 and Planck Collaborations], “A Joint Analysis of BICEP2/Keck 
Array and Planck Data,” Submitted to: Phys.Rev.Lett. [arXiv:1502.00612 [astro-ph.GO]]. 

[8] A. H. Guth, “The Inflationary Universe: A Possible Solution to the Horizon and Flatness 
Problems,” Phys. Rev. D 23, 347 (1981). A. D. Linde, “A New Inflationary Universe Scenario: 
A Possible Solution of the Horizon, Flatness, Homogeneity, Isotropy and Primordial Monopole 
Problems,” Phys. Lett. B 108, 389 (1982). A. Albrecht and P. J. Steinhardt, “Cosmology for 
Grand Unified Theories with Radiatively Induced Symmetry Breaking,” Phys. Rev. Lett. 48, 
1220 (1982). 

[9] A. D. Linde, “Chaotic Inflation,” Phys. Lett. B 129, 177 (1983). 

[10] R. Allahverdi, R. Brandenberger, F. Y. Cyr-Racine and A. Mazumdar, “Reheating in 
Inflationary Cosmology: Theory and Applications,” Ann. Rev. Nucl. Part. Sci. 60, 27 (2010) 
[arXiv: 1001.2600 [hep-th]]. 

[11] J. Martin and C. Ringeval, JCAP 0608, 009 (2006) [astro-ph/0605367]. 

[12] J. Martin and C. Ringeval, “First CMB Constraints on the Inflationary Reheating 
Temperature,” Phys. Rev. D 82 , 023511 (2010) [arXiv: 1004.5525 [astro-ph.GO]]. 

[13] J. Martin, C. Ringeval and V. Vennin, “Encyclopdia Inflationaris,” Phys. Dark Univ. (2014) 
[arXiv:1303.3787 [astro-ph.GO]]. 

[14] L. Dai, M. Kamionkowski and J. Wang, “Reheating constraints to inflationary models,” Phys. 
Rev. Lett. 113, 041302 (2014) [arXiv: 1404.6704 [astro-ph.GO]]. 

[15] M. Drewes and J. U. Kang, “The Kinematics of Cosmic Reheating,” Nucl. Phys. B 875, 315 
(2013) [Nucl. Phys. B 888, 284 (2014)] [arXiv:1305.0267 [hep-ph]]. M. Drewes, “On finite 
density effects on cosmic reheating and moduli decay and implications for Dark Matter 
production,” JCAP 1411, no. 11, 020 (2014) [arXiv: 1406.6243 [hep-ph]]. 

[16] V. A. Belinsky, 1. M. Khalatnikov, L. P. Grishchuk and Y. B. Zeldovich, “Inflationary Stages 
In Cosmological Models With A Scalar Field,” Phys. Lett. B 155, 232 (1985). 

[17] P. Creminelli, D. Lpez Nacir, M. Simonovi, G. Trevisan and M. Zaldarriaga, “ffi or Not (jfi : 
Testing the Simplest Inflationary Potential,” Phys. Rev. Lett. 112, 241303 (2014) 

[arXiv: 1404.1065 [astro-ph.GO]]. P. Creminelli, D. L. Nacir, M. Simonovi, G. Trevisan and 
M. Zaldarriaga, “(jfi Inflation at its Endpoint,” arXiv:1405.6264 [astro-ph.GO]. 

[18] L. McAllister, E. Silverstein and A. Westphal, “Gravity Waves and Linear Inflation from 
Axion Monodromy,” Phys. Rev. D 82, 046003 (2010) [arXiv:0808.0706 [hep-th]]. 

[19] M. S. Turner, “Coherent Scalar Field Oscillations in an Expanding Universe,” Phys. Rev. D 28 
(1983) 1243. 


- 15 - 



[20] L. F. Abbott, E. Farhi and M. B. Wise, “Particle Production in the New Inflationary 
Cosmology,” Phys. Lett. B 117, 29 (1982). 

[21] A. D. Dolgov and A. D. Linde, “Baryon Asymmetry in Inflationary Universe,” Phys. Lett. B 
116, 329 (1982). 

[22] A. Albrecht, P. J. Steinhardt, M. S. Turner and F. Wilczek, “Reheating an Inflationary 
Universe,” Phys. Rev. Lett. 48, 1437 (1982). 

[23] D. I. Podolsky, G. N. Felder, L. Kofman and M. Peloso, “Equation of state and beginning of 
thermalization after preheating,” Phys. Rev. D 73, 023501 (2006) [hep-ph/0507096]. 

J. F. Dufaux, G. N. Felder, L. Kofman, M. Peloso and D. Podolsky, “Preheating with trilinear 
interactions: Tachyonic resonance,” JCAP 0607, 006 (2006) [hep-ph/0602144]. 

[24] V. F. Mukhanov and G. V. Chibisov, “Quantum Fluctuation and Nonsingular Universe. (In 
Russian),” JETP Lett. 33, 532 (1981) [Pisma Zh. Eksp. Teor. Fiz. 33, 549 (1981)]. 

S. W. Hawking, “The Development of Irregularities in a Single Bubble Inflationary Universe,” 
Phys. Lett. B 115, 295 (1982). A. A. Starobinsky, “Dynamics of Phase Transition in the New 
Inflationary Universe Scenario and Generation of Perturbations,” Phys. Lett. B 117, 175 
(1982). A. H. Guth and S. Y. Pi, “Fluctuations in the New Inflationary Universe,” Phys. Rev. 
Lett. 49, 1110 (1982). J. M. Bardeen, P. J. Steinhardt and M. S. Turner, “Spontaneous 
Creation of Almost Scale - Free Density Perturbations in an Inflationary Universe, ” Phys. Rev. 
D 28, 679 (1983). 

[25] P. A. R. Ade et al. [Planck Gollaboration], “Planck 2015 results. XIII. Cosmological 
parameters,” arXiv:1502.01589 [astro-ph.GO]. 

[26] J. B. Munoz and M. Kamionkowski, “Equation-of-State Parameter for Reheating,” Phys. Rev. 
D 91, no. 4, 043521 (2015) [arXiv:1412.0656 [astro-ph.GO]]. J. Martin, G. Ringeval and 

V. Vennin, “Observing Inflationary Reheating,” Phys. Rev. Lett. 114, no. 8, 081303 (2015) 
[arXiv:1410.7958 [astro-ph.GO]]. J. Martin, “The Observational Status of Cosmic Inflation after 
Planck,” arXiv: 1502.05733 [astro-ph.GO]. J. L. Gook, E. Dimastrogiovanni, D. A. Easson and 
L. M. Krauss, “Reheating predictions in single field inflation,” arXiv:1502.04673 [astro-ph.GO]. 

[27] V. Domcke and J. Heisig, “Constraints on the reheating temperature from sizable tensor 
modes,” arXiv: 1504.00345 [astro-ph.GO]. 

[28] J. Beringer et al. (Particle Data Group), Phys. Rev. D86, 010001 (2012) 


- 16 - 



